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Abstract 

We investigate the asymptotic behavior as time goes to infinity of Hawkes processes whose 
regression kernel has norm close to one and power law tail of the form with 

a G (0,1). We in particular prove that when a G (1/2,1), after suitable rescaling, their 
law converges to that of a kind of integrated fractional Cox-Ingersoll-Ross process, with 
associated Hurst parameter H = a — 1/2. This result is in contrast to the case of a 
regression kernel with light tail, where a classical Brownian CIR process is obtained at 
the limit. Interestingly, it shows that persistence properties in the point process can lead 
to an irregular behavior of the limiting process. This theoretical result enables us to give 
an agent-based foundation to some recent findings about the rough nature of volatility in 
financial markets. 

Keywords: Hawkes processes, limit theorems, nearly unstable processes, heavy tail, frac¬ 
tional stochastic equation, fractional Cox-Ingersoll-Ross process, volatility, long memory. 


1 Introduction 


A Hawkes process {Nt)t>o is a self-exciting point process whose intensity at time t, denoted 
by Ai, is of the form 


At = p -b X 

0<Ji<t 


p + 



(j){t - s)dNs, 


where p is a positive real number, </) a non-negative measurable function and the Ji are the 
points of the process before time t (see Section 2 for a more formal definition). These pro¬ 
cesses have been introduced in the early seventies by Hawkes, see [21, 22, 23], in the purpose 
of modeling earthquakes and their aftershocks, see [1] for such application. In the last years, 
the probabilistic and statistical analysis of Hawkes processes has known several interesting 
developments, driven by the recent use of Hawkes processes in various applied fields such as 
neurosciences [13, 33, 34, 36], sociology [9, 28, 41], criminology [30, 31], genome analysis [37] 
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and mostly finance [2, 4, 7, 10, 12, 14, 15]. 


Among the probabilistic questions raised by Hawkes processes, particular attention has been 
devoted to the study of their long term scaling limits. More precisely, one wishes to understand 
the behavior as T tends to infinity of the process 


ariNtT), t G [0,1], 

where ot is a suitable normalizing factor. In [5], it is shown that under the condition 


1 — 


r+oo 

/ 4>is)d 
Jo 


s < 1, 


the asymptotic behavior of a Hawkes process is quite similar to that of a Poisson process. 
Indeed, as T tends to infinity, 


TG'r^i'TlI .n 

sup |— -E[—J| ^ 0, 
4G[0,1] ^ ^ 


in probability and 


NtT r J^tT 

— -E — 

rp L rp 


tg[0,l] 




in law for the Skorohod topolgy, with a an explicit constant and [Wt) a Brownian motion. 
This result has been extended in [42] to the case of non-linear Hawkes processes. 


The condition < 1 is essential in order to obtain the preceding result. It is actually very 
similar to the assumption ]/)] < 1 one makes on the autoregressive coefficient p when working 
with a discrete time stationary AR(1) process. In particular, when starting the Hawkes pro¬ 
cess at t = —oo, the assumption < 1 is required in order to get a stationary intensity 
with finite hrst moment. Also, as for AR(1) processes, under this condition, Hawkes processes 
only exhibit weak dependence properties. Consequently, their asymptotic behavior is in that 
case no surprise, close to that of a Poisson process. Hence this condition is called stability 
condition. 


In [26], the authors investigate the scaling limit of Hawkes processes when the stability con¬ 
dition is almost violated. This means they consider a sequence of Hawkes processes satisfying 
the stability condition, but for which the kernel (f) = cfP" also depends on the observation scale 
T, such that tends to 1 as T goes to infinity. Such a sequence is called sequence of 

nearly unstable Hawkes processes. 

Beyond its obvious mathematical interest, considering the case of nearly unstable Hawkes 
processes is motivated by empirical studies on financial data. Indeed, it has become quite 
standard to model the clustered nature of order flows on financial markets by means of Hawkes 
processes. However, one systematically estimates norms for the regression kernels which 
are smaller but very close to 1, see [16, 17, 19, 27]. Interestingly, this empirical stylized fact 
that Hawkes processes have to be nearly unstable to fit the data has a very natural financial 
interpretation, namely the high degree of endogeneity of modern markets due to high fre¬ 
quency trading. This signifies that a large proportion of orders is just endogenously triggered 
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by other orders, see [17, 19, 26] for more details. In this framework, it is proved in [26] that 
the limiting law of a sequence of nearly unstable Hawkes processes is that of an integrated 
Cox-Ingersoll-Ross process (CIR process for short). Hence, compared to the case where the 
stability condition is in force, the asymptotic behavior at first order is no longer deterministic, 
see also [42] for the case where ||())||i is exactly equal to one and other interesting developments. 
Note that this CIR scaling limit seems to be very consistent with hnancial practice. Indeed, 
it is widely acknowledged that there exists a linear relationship between the cumulated order 
flow and the integrated squared volatility, see for example [40], and CIR processes are very 
classical models for the squared volatility. 


Nevertheless, the CIR limit in law of nearly unstable Hawkes processes discussed above is 
obtained under the crucial assumption 

r+oo 

/ s(l){s)ds < -boo. 

Jo 

It is therefore quite natural to try to extend the results of [26] to the case of nearly unstable 
heavy tailed Hawkes processes, for which this condition is no longer satisfied. Hence we 
consider in this paper the situation where 


(j){x) 


K 

rsj - 

x—>-+oo 


where a G (0,1) and K \s a, positive constant. This setting is actually much more in agreement 
with hnancial data, where one not only hnds that the function (j) has an norm close to one, 
but also that it has a power law tail, see [6, 19]. This heavy tail is quite easy to interpret 
in practice too: it is related to the persistence of the signed order how (the series of -|-1, —1 
where -|-1 represents a buy order and —la sell order). Indeed, the long memory property of 
this process is well established and is due to the so-called order splitting phenomenon: most 
orders are actually part of large orders (called metaorders), which are split in smaller orders 
so that prohibitive execution costs can be avoided. 


Our main result is that for a G (1/2,1), after proper rescaling, the law of a nearly unstable 
heavy tailed Hawkes process converges to that of a process which can be interpreted as an 
integrated fractional diffusion. Loosely speaking, this limiting distribution can be viewed as 
the integral of a fractional version of the CIR process, where a fractional Brownian motion 
replaces the ordinary Brownian motion. This result is quite remarkable from a probabilis¬ 
tic point of view. Indeed, assuming fat tail leads to a limit which is not an integrated 
semi-martingale. This is in strong contrast to all other scaling limits obtained for Hawkes 
processes. Technically, this heavy tail case is of course more subtle than that investigated in 
[26] where semi-martingale theorems are used in a quite direct manner. Moreover, Gaussian 
methods are not easy to apply in our context since the limit is not a simple Gaussian func¬ 
tional, although it somehow involves a fractional Brownian motion. 


The perhaps most surprising phenomenon obtained in our result is the value of the Hurst 
parameter H of the (sort of) fractional Brownian motion appearing in the limit. Indeed, 
fat tail meaning persistence, one would expect getting also persistence in the limit and so 
H > 1/2. This is actually the contrary: an aggregation phenomenon occurs in the heavy tail 
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case, leading to a very irregular process in the limit, its derivative behaving as a fractional 
Brownian motion with Hurst parameter H < 1/2. Coming back to financial applications, this 
means that in practice, the volatility process should be very irregular, which is perfectly in 
line with the recent empirical measures of the volatility smoothness obtained in [18]. There¬ 
fore, our theoretical result shows quite clearly that the rough behavior of the volatility can 
be explained by the high degree of endogeneity of financial markets together with the order 
splitting phenomenon. This is to our knowledge the first agent-based explanation for the very 
rough nature of the volatility. 

The paper is organized as follows. We first give our assumptions together with some intuitions 
about the limiting behavior of our processes in Section 2. Section 3 contains our main theorems 
whose proofs can be found in Section 4. Finally, some technical results are relegated to an 
appendix. 


2 Assumptions and intuitions for the results 

We describe in this section our asymptotic framework together with intuitions about our main 
results which are given in Section 3. 


We consider a sequence of point processes indexed by T^. For a given T, {Nf) 

satisfies Nq = 0 and the process is observed on the time interval [0, Tj. Our asymptotic 
setting is that the observation scale T goes to infinity. The intensity process {Xj) is defined 
for t > 0 by 

>H = - s)dNj, 

Jo 

where is a sequence of positive real numbers and the are non-negative measurable 
functions on M'*' which satisfies < -|-oo. For a given T, the process {Nf) is defined 

on a probability space equipped with the filtration (-T7')tG[o,T]) where is 

the cj-algebra generated by {Nj)s<t- Moreover we assume that for any 0 < a < b < T and 

niNl-Nj)U]=n[ AjUds], 

J a 


which sets as the intensity of N'^. In particular, if we denote by {Jn)n>i the jump times 
of {NJ'), the process 




X'^ds 


is a martingale and the law of N'^ is characterized by A^. From [24], such a construction can 
be done and the process is called a Hawkes process. 


Let us now give more specific assumptions on the functions (jx^. 
Assumption 1. For t G ]!“*■, 

= aT4>{t), 

^Of course by T we implicitly means Tn with n G N tending to infinity. 
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where (aT)T>o is a sequence of positive numbers converging to 1 such that for all T, ot < ^ 
and 4> is a non-negative measurable function such that ||(/>||i = 1. Furthermore, 


lim ax“(l 

X—>- + 00 


- F{x)) = K, 


for some a G (0,1) and some positive constant K, with 


F{x) = f (l){s)ds. 

Jo 

Recall that in [26], it is assumed that 



t(l){t)dt 


< +00 


( 1 ) 


and this condition leads to a CIR-type limit. Considering Assumption 1 instead of (1) will 
induce a completely different scaling behavior for the sequence of nearly unstable Hawkes 
processes. Nevertheless, in this framework, we still have almost surely no explosion^; 


lim jJ = + 00 . 

rx^+oo 


Remark that we do not work in the stationary setting since our processes start at time t = 0 
and not at t = —oo. 


Let denote the martingale process associated to that is, for f > 0, 

Mf = Nf - f \Us. 

Jo 

We also set as the function defined on by 

OO 

^^(t) = ( 2 ) 

k=l 

where {(jJ^)*^ = and for k >2, denotes the convolution product of with 

the function (jp". Note that pP"{t) is well defined since < 1- This function plays an 

important role in the study of Hawkes processes, see [3]. In particular, it is proved in [26] 
that the intensity process, rescaled on [0,1], can be rewritten 

ftT ftT 

Xfrp = fp "/ p'^{Tt — s)pP’ds / p'^{Tt — s)dMj. 

Jo Jo 

In term of scaling in space, a natural multiplicative factor is (I — aT)/pP'- Indeed, in the 
stationary case, the expectation of Xj is pP /{t — lli;i)^lli). Thus, the order of magnitude of 
the intensity is pPp — otP^- This is why we define 

riT _ (1 “ ^t) 

^In fact, for a Hawkes process, the no explosion property can be obtained under weaker conditions, for 
example <j){s)ds < oo for any t > 0, see [5]. 
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Then we easily get 


Cj' = {l-aT)+[ T{1 - aT)ip'^{Ts)ds + [ ip'^{T{t-s)) 

Jo \j d- Jo 


CTdBj 


( 3 ) 


with 



From (3), we see that the asymptotic behavior of the intensity is closely related to that of 
X 'ifF{Tx). To analyze the limiting behavior of this function, let us remark that for x > 0, 


p^{x) = T 




is the density of the random variable 



i=l 


( 4 ) 


where the (Xj) are iid random variables with density cj) and is a geometric random variable 
with parameter 1 — The Laplace transform of the random variable denoted by p^, 
satisfies for z >0 


f{z)=E[e-^‘^^] = ^(1 -aT)(oT)^“^IE[e-^^^-i^*] 


k=l 


£(l-aT)(aT)‘-‘{?(^))‘ 


</>( 


.=1 - 

where (p denotes the Laplace of cp. We now need to compute an expansion for (p{z). Using 
integration by parts, we get 

^ p-\-OQ ^+00 

(p{z) = z / e~^^F{t)dt = 1 — z / — F{t))dt. 

Jo Jo 

Then using Assumption 1 together with Karamata Tauberian theorem (see for example The¬ 
orem 17.6 in [8]), we get 

^{z) = 1 - + o(z“), 

a 

with T the gamma function. Set 6 = and vt = — a-r). As T goes to infinity, 

(z) is thus equivalent to 

(5) 

Vt + 2“ 

The function whose Laplace transform is equal to this last quantity is given by 

VTX'^~^Ea,a{-VTX°‘), 


> 0, P[/^ = fc] = (1 - aFiarY 
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with the (a, /3) Mittag-Leffler function, that is 


Ea^piz) — 


^-r(an + /3)’ 


see [20]. Putting this together with (3) and (4), we can expect (for a > 1/2) 

Cj ~ Ur / s'"~^Ea,a{-VTS°')ds + -fTVT [ {t - Ea,a{-VT{t - s)"") dB^, 
Jo Jo 


with 


1 


7t = 


\//i^r(l - ar) 

The process can be shown to converge to a Brownian motion B. Thus, denoting by Uoo 
and 7oo the limits of vt and 7 ^, passing (non rigorously) to the limit, we obtain (for a > 1/2) 




[ s“ ^Ea,a{-VooS°')ds+ 'JooVoo [ (^ “ s)“ ^Ea,a{-Voo{t - s)'^)^/C^dBs. (6) 

Jo Jo 


From (6), we see that in order to get a non-deterministic asymptotic behavior for Cj, we 
need that both Voo and 7 cxd are positive constants or Voo is equal to zero and 'jooVoo is positive. 
However, in the last situation, the expectation of (C/^)^ would be of order Such cases 

where the variance is increasing with power law rate are incompatible with the approximate 
stationarity property we want to keep for our model and its limit. Therefore, only one regime 
seems to be natural and this leads us to the following assumption. 

Assumption 2. There are two positive constants A and p* such that 


and 


lim r"(l -ut) = A<5. 

T^+oo 


lim ri-“uT = 
T—>-+oo 


In particular. Assumption 2 implies that vt converges to A and therefore the sequence of 
random variables (J^) converges in law towards the random variable whose density on is 
given by 

Aa:"-^F;a,„(-Ax“). 

Beyond giving the suitable asymptotic regimes for ax and ht, the heuristic derivation leading 
to (6) provides us an expression for the limiting law of the rescaled intensities of our sequence 
of nearly unstable heavy tailed Hawkes processes. In (6), this law appears under the form of 
a non-classic stochastic integral equation. Indeed, it is of Volterra-type and is therefore (a 
priori) neither a diffusion nor a semi-martingale. Furthermore, the main term of the Volterra 
kernel x°‘~^ exhibits a singularity at point 0, of the same kind as that of the fractional 
Brownian motion (B^^) when expressed under the form: 


= 


I 


F(F + I/2) Vo 


{ [\t-sf-^/‘^dWs+ f {t-s)^-^/^-{-s)^-^/^dWs), 

J —oo 


(7) 


7 







with {Wt) a Brownian motion, see [29]. 

The preceding computations suggest a possible approach to derive the limiting behavior of our 
sequence of Hawkes processes: studying the intensity of the processes. Indeed the intensities 
can be rewritten under the form of stochastic integral equations as (3). Consequently, one 
can try to pass to the limit in the coefficients of the equation to obtain the limiting law, as 
we (non rigorously) did to get (6). This is exactly the approach used in [26]. However, in this 
more intricate heavy tail case, it seems very hard to use. In particular, the sequence {Cf) 
is typically not tight. Thus, instead of considering the intensities, we directly work on the 
Hawkes processes themselves, more in the spirit of [42]. 


3 Main results 

We rigorously state in this section our theorems on the limiting behavior of nearly unstable 
heavy tailed Hawkes processes. We start with some technical results about the function 
appearing as the inverse Laplace transform of (5) in Section 2. 


3.1 The function 

As shown by the derivations in the previous section, the function 

/“A(x) = 


plays a crucial role in our analysis. We give here some elements about the regularity of 
this function which will be useful in the sequel. We denote by /“/ and D^f the fractional 
integration and derivation operators, which are defined for a suitable measurable function / 
by 

fit) 




and 


D'^fix) = 


{x — 

d r fit) 


dt 


-dt. 


r(l — a) dx Jq {x — t)'^ 

The following lemma is a direct consequence of the definition of and Section 11 in [20]. 
Proposition 3.1. The function f°‘’^ is C°° on (0,1] and 


Furthermore, 
For V < a, f 


r\x) 


A 


x—>-0+ r(Q;) 


X 


a —1 


(/“’')'(x) 




a:->-0+ r(Q:) 

has Holder regularity a on (0,1]. 

is V fractionally differentiable and 

D-^rix) = 


8 








Therefore, 


and 


D’^rix) 
{D’^rrix) ~ 


A 


x^o+ r(a - v) 

A(a — 1 — v) 1 

x->-o+ r(a — v) 


For v' > 0, /“ is v' fractionally integrable and 


1 


rix) = A^3^i?„,„+,.(-Ax“). 


X 


Therefore, 


and for a + z/' / 1, 


r r (x) 


A 


1 


(r'r)'(x) 


;^0+ r(a + u') x^~°‘~’^' 

A(q! — 1 + ly') 1 


>0+ T{a + v') 

Proposition 3.1 will be a key tool in the proofs of the main results. 


3.2 The limiting behavior of nearly unstable heavy tailed Hawkes processes 

Let us first give some notations. We consider for t G [0,1] the renormalized Hawkes process 


Xf = 


1 — ax 


T^Ijl*6 


m 




and its associated integrated intensity 


Af = 


1 — ax 




X^ds. 


As explained in Section 2, the space renormalization is chosen so that the processes have an 
expectation of order one. We also introduce the martingale defined on [0,1] by 


zl 


1 — ax 


ixl 



We are now ready to give our results about the convergence in distribution of for 

the Skorohod topology. 

Proposition 3.2. Under Assumptions 1 and 2, the sequence {Z'^,X'^) is tight. Furthermore, 
if {Z,X) is a limit point of {Z'^,X'^), then Z is a continuous martingale and [Z,Z] = X. 

Now let (Z,X) be a couple of processes defined on some probability space (H,Al,P) with 
law being one of the possible limit points of the sequence of distributions associated to the 
sequence {Z'^,X'^). From Proposition 3.2, we are able to obtain the following theorem which 
is one of our main results. 
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Theorem 3.1. There exists a Brownian motion B on (r2,.4., P) (up to extension of the space) 
such that for t G [0,1], Zt = Bxt and for any e > 0, X is continuous with Holder regularity 
(1 A 2a) — e on [0,1] and satisfies 

Xt = sr’\t - s)ds + r’\t - s)BxJs. (8) 

Hence the limiting process in Theorem 3.1 has a quite original form, which can actually be 
interpreted more easily by looking at its derivative (when it exists). 

3.3 The limiting volatility process 

As explained in the introduction, when it exists, the derivative of the limiting process X in 
Theorem 3.1 can be interpreted as a volatility function. Actually, if the tail of the function (j) 
is not too heavy, X is indeed differentiable. Let us write 

= [ r’^{s)ds. 

Jo 

The following result holds. 

Theorem 3.2. Let (Xt) be a process satisfying (8) fort G [0,1] and assume a > 1/2. Then 
X is differentiable on [0,1] and its derivative Y satisfies 

Yt = F^’\t) + r\t - s)^sdBl (9) 

with B^ a Brownian motion. Furthermore, for any e > 0, Y has Holder regularity a —1/2 —e. 

3.4 Discussion 

We now comment the results given in Theorem 3.1 and Theorem 3.2. 

• The singularity at zero of the function /“A appearing in our two theorems is of order 
Making an analogy with the Volterra representation of the fractional Brownian motion (7), 
this corresponds to a Hurst parameter H equal to a — 1/2. Thus, in the case a > 1/2 where 
our volatility process is well defined, because of the square root term in front of the Brownian 
motion, we can somehow interpret (9) as a fractional CIR process with Hurst parameter equal 
to a — 1/2. This dynamics leads to a very rough process, with Holder regularity close to zero 
when a is close to 1/2. As mentioned in the introduction, this is perfectly consistent with 
recent empirical measures of the volatility smoothness on financial data, see [18]. 

• A practical consequence of the preceding point is the following: When observing on a time 
interval of order 

1 

(i-ll-MliA" 

a Hawkes process with kernel (j) with norm close to one and power law tail with index 
1 + ce, then after rescaling, a fractional-like behavior is obtained. 
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• Theorem 3.2 relates the smoothness of the volatility process to the tail parameter a. 
This is particularly interesting for financial applications. Indeed, the parameter a is usu¬ 
ally considered very hard to measure. Our theorem provides an approach where it can be 
obtained relying on the smoothness of the volatility, which is much easier to estimate, see [18]. 

• The irregular volatility appearing at the limit arises because our Hawkes processes are 
nearly unstable with heavy tailed kernels. As explained in the introduction, in financial 
terms, it means that the rough behavior of the volatility can be explained by the high de¬ 
gree of endogeneity of modern markets combined with the persistent nature of the order flows. 

• A natural question is that of the uniqueness of the solution of Equation (9). Because 
of the singular Volterra kernel and of the square root term, it is probably quite difficult to 
answer. Actually, in the very recent paper [32], using SPDE techniques, the authors show 
weak uniqueness for the solutions of an equation which is quite similar to (9). However, they 
use the specific form of their equation and their approach cannot be adapted to our case in 
an obvious way. 

• Note that Hawkes processes with norm exactly equal to one have been introduced in 
[11]. In this work, the authors show that in order to get a stationary intensity, the parameter 
[1 must be equal to zero and the regression kernel has to be heavy tailed. Several additional 
results for the non-stationary heavy tailed case (||0||i = 1 and fJ. > 0) can be found in [42]. 

• Compared to the approach in [26], it is important to remark that our volatility process 
is simply the derivative of the limit of the sequence of nearly unstable heavy tailed Hawkes 
processes. Contrary to what is done in [26], we do not provide any result about the convergence 
of the sequence of intensities of the Hawkes processes. In particular, the sequence of intensities 
is not shown to converge towards the volatility^. Note also that our assumptions are slightly 
weaker than those in [26]. In particular, we do not require the function cf) to be bounded. 
Again, this is relevant for financial applications where becomes typically very large as t 
tends to zero, see [6]. 

4 Proofs 

We give in this section the proofs of Proposition 3.2, Theorem 3.1 and Theorem 3.2. In the 
sequel, c denotes a positive constant which may vary from line to line (and even within the 
same line if no ambiguity). 

4.1 Proof of Proposition 3.2 

We show here the tightness of We start with the following lemma. 

Lemma 4.1. The sequences and A^ are C-tight. 

Proof. Erom [5], we get that the expectation of the Hawkes process Nf satisfies 

lE[iVf’] = + PT f - s)sds < + IIV^'^lli)- 

Jo 

Actually it can be shown that for some reasonable functions (j), the sequence of intensities does not converge, 
at least in the Skorohod topology. 
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Therefore, since 


c 


we get 



1 ’ 
1 — ax 


E[Xf] = E[Af ] < c. 


The tightnesses of and then follow, using the fact that both processes are increasing. 


Moreover, since (1 — aT)/T°‘ tends to zero, the maximum jump size of X'^ and A^ (which is 
continuous) goes to zero as T tends to infinity. From Proposition VI-3.26 in [25], this implies 
the C-tightness of X'^ and A"^. □ 

We now give the proof of Proposition 3.2. It is easy to get that the angle bracket of is 
A^. Prom Lemma 4.1, it is C-tight. Thus, from Theorem VI-4.13 in [25], the sequence {Z'^) 
is tight. Finally, marginal tightnesses imply the joint tightness of (Z^,X^). 

Let us now consider a subsequence converging towards a process that we denote 

by {Z,X). Using Proposition VI-6.26 in [25] together with the fact that the bracket of 
is we get that X = [Z,Z]. 

Since goes to 0, the maximum jump size of tends to zero. Therefore, is 

C-tight and so the limit Z is continuous. It remains to show that Z is a martingale. Using 
Corollary IX. 1.19 in [25], Z is a local martingale. Moreover, the expectation of its bracket 
being finite, it is a martingale. 


4.2 Proof of Equation (8) in Theorem 3.1 

We start with the following lemma which shows that we can somehow work with A^ rather 
than with X'^. 


Lemma 4.2. The sequence of martingales X^ — h?" tends to zero in probability, uniformly 
on [0,1]. 


Proof. We have 


X/ - A/ = 


s:_l 


T°‘n*S- 


Applying Doob’s inequality to the martingale M^, we get 


E[ sup {{Xj’ - Aff}] < c( ^ ^^ )^E[(Mf)^]. 

te[o,i] J 


Then, the bracket of being X^, we deduce 


E[ sup {{Xj' 

tG[0,l] 


Af)"}] < c (—J 


'J'2c 


which ends the proof. 

We now state a lemma which will be useful in the proof of Equation (8). 


□ 
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Lemma 4.3. The sequence of measures with density p^{x) defined by Equation (4) converges 
weakly towards the measure with density Ax“). In particular, over [0,1], 


converges uniformly towards 



Proof. The proof of this result is obtained by showing that the Laplace transform of the 
measure with density p^{x) converges towards the Laplace transform of the measure with 
density Xx°^~^Ea,a{—^x°^). This has already been done in Section 2. □ 


We now give the proof of Equation (8). Let us consider a converging subsequence , X'^^) 
and write {Z, X) its limit. Abusing notation slightly, we write {Z"^, X'^) instead of (Z'^”, X'^"). 
Using Skorokhod’s representation theorem, there exists a probability space on which one can 
define copies in law of the {Z'^, X^) converging almost surely for the Skorohod topology to a 
random variable with the same law as (Z,X). We now work with this sequence of variables 
converging almost surely and their limit. The processes Z and X being continuous, we have 


sup \X'[ — Xt\^ 0, sup \Z'[ — Zt\ ^ 0. (10) 

tG[0,l] tG[0,l] 


Let us now rewrite the cumulated intensity. For all t > 0, we have 

rt rt rs pt 

/ X^ds = tfiT+ / — X^du)ds + / {t — s)Mj'ds. 

Jo Jo Jo Jo 

Then, using that * cjfi" = — cjfi", where is defined in Equation (2), remark that 

pt pS pt pt 

/ fi'^{t — s) / (p^{s — r)Mj'drds = / / lr<sfi^{t — s)^^{s — r)dsMjdr 

Jo Jo Jo Jo 


n t—r 

[t — r — s)4J" {s)dsM^ dr 

= f * (fp"{t — r)Mj'dr 
Jo 

= f fP"{t — r)Mjdr — f cjp{t — r)Mjdr. 

Jo Jo 

This together with Lemma 3 in [5] yields 

f Xjds = tfiT + f fP {t — s)syLTds + f fP {t — s)mJ ds. 

Jo Jo Jo 

Therefore, replacing t by Tt, multiplying by (1 — ot) /{T°^ fi*5~p, and writing 

fiT 


ut = 
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we get A'^(t) = Ti + T 2 + T 3 , with 


Ti = (1 — aT)tuT, 

T2 = T{l — aT)uT / (T{t — s))sds, 
Jo 

rp^ _ pl—af2 . / (^ ) 




ds. 


Since ut converges to 1, we get that Ti goes to zero. For T 2 , note that integrating by parts, 
we have 


T2 = 


cltut / — s)sds = utut / F'^{t — s)ds. 

Jo Jo 


Using Lemma 4.3 and integrating by parts again, we obtain that T 2 tends uniformly to 


f F^’^{t-s)ds= [ - s)sds. 

Jo Jo 


We now turn to T 3 . Remark that 


T3 = 


Qt 


- aT)f^*d- 


=f Jo P^{t-s)Zjds 


and recall that 


zl = 




1 — Qt 


-{Xt -Ai). 


Thus, using that is piecewise constant, applying integration by parts, we get (pathwise) 


and in the same way 


Then, 


[ p^{t-s)Zfds= [ F^{t-s)dZf 
Jo Jo 

[ s)Z^ds= [ F^’^{t- s)dzj. 

Jo Jo 


e\{ [ {F^’^{t 
L Jo 

which tends to zero thanks to Lemma 4.3. Furthermore, using (10), we get that 


- s) - F^{t - s))dZff <c -s)- F^{t - s)fds, 

-I 40 


f\{t 

Jo 


- - sr)\Zs - Zl \ds 


also tends to zero. Consequently, we finally obtain that for any t, T 3 converges to 


1 


y/pF\ 


\{t - - sY)Zsds. 


Since Z is a continuous martingale, the fact that Zt = Bxt is a consequence of the Dambis- 
Dubin-Schwarz theorem, see for example Theorem V-1.6 in [35]. 
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4.3 Proof of the Holder property for X in Theorem 3.1 

We start with the following lemma. 

Lemma 4.4. Let B he a Brownian motion and X a solution of (8) associated to B. Let 
H in (0,1). If X has Holder regularity H on [0,1], then for any e > 0, X has also Holder 
regularity {{a + H/2) A 1) — e on [0,1]. 

Proof. Let e > 0 and Zt = Bxf The function 

t ^ [ sf^’^{t — s)ds 

Jo 

being C^, it is enough to show that 

[ r’^{t-s)Zsds 
Jo 

has Holder regularity ((a + H/2) A 1) — e. Since for any e' > 0, Z has Holder regularity 
{H/2 — e'), by Proposition A.l, it is {H/2 — e) fractionally differentiable and is 

continuous. Using the fact that is fractionally integrable, from Corollary A.l, we get 

[ r’^{t - s)Zsds = [ l^l’^-^f^'^{t-s)D^I‘^-^Z,ds. 

Jo Jo 

Finally, the properties of /W 2 -£y‘a,A Proposition 3.1 together with Proposition A.3 

give the result. □ 

We now show that for H be a Brownian motion and X a solution of ( 8 ) associated to B, 
then, for any e > 0, almost surely, the process X has Holder regularity (1 A 2a) — e on [0,1]. 

Let M be the supremum of the Holder exponents of X. From Proposition 3.1 together with 
Proposition A.3, we get that M > a. 

Let us now assume that M < (1 A 2a). Then we can find some H < M and some e > 0 such 
that 

M < ((a + H/2) A 1) - e. 

Thus, since X has Holder regularity H, Lemma 4.4 implies that X has also Holder regularity 

((o + H/2) A l) — e, 

which is a contradiction. Therefore M > (1 A 2a), which ends the proof. 

4.4 Proof of Theorem 3.2 

First remark that thanks to the Holder property of the Brownian motion together with that 
of the process X, we immediately deduce the following lemma. 

Lemma 4.5. Let B be a Brownian motion, X a solution of (8) associated to B and Zt = Bxf 
Then, for any e > 0, almost surely, the process Z has Holder regularity {1/2 f\a) — e on [0,1]. 
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We now give the proof of Theorem 3.2. Using Proposition 3.1, Lemma 4.5 and Corollary A.2 
for any v G (0, a), we can rewrite Equation (8) as 


= r sn\t - s)ds + ^1= f - s)FZsds. 

Jo V d- ^ Jo 

Moreover, taking > 1/2, since Z is 1 — fractionally differentiable, we get 

rZs = [ D^-'^Zudu. 

Jo 


Thus, using Fubini’s theorem, we obtain 



Hence, we get 



Zududs 
s)D^~‘' Zudsdu 
Zudsdu 
Zududs. 


with 

W = r - u)D^-'^Zudu. 

V/^*A Jq 

From Proposition 3.1 together with Proposition A.3, we have that Y has Holder regularity 
{a — v). Thus, taking v close enough to 1/2, we get that for any e > 0, T has Holder regularity 
(a — 1/2 — e). This implies that X is differentiable with derivative Y. 


Now, since Z is a continuous martingale with bracket X and because v > 1/2, we can use 
the stochastic Fubini theorem, see for example [39], to obtain 


Therefore, 




1 

d 

r 

r(i. 

) ds 

Jo 

1 

d 

r 

r(z. 

) ds 

Jo 

1 

d 

r 

r(z. 

) ds 

Jo 


1 

d 


Z,, 


-dv 


\l-u 


dZudv 

dvdZ,, 


r(l^ + 1) (is Jq 


[\s-u) 

Jo 


- uYdZ„,. 


Yt = + 


y/jYX 


D 


'nt-s) [ {s-ufdZuds. 

T{u + 1) ds Jo 
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Using Fubini’s theorem twice and the fact that / * {g') = {f * g)', we derive 

= F»’\t) + r{t - u)dZu. 

Moreover, using Theorem V-3.8 of [35], there exists a Brownian motion such that 

f ./YsdBl 

Jo 

So, consider now the process (Yj) dehned by 


Yt = ^ r{t - u)^/Y,dBl. 


Going backward in the previous computations for Yt and Zg, we remark that 


Yt = F^’\t) + 


1 


1 


d 




[ {s - uY s/YidBl^ds 
Jo 


and 


1 d 
T{v + 1) ds 


f‘(s - urVY^dBl = j-Y 

1 d 


r(i/) ds Jo (■s — 

1 d r z, 


V^dBi)dv 


r(i/) ds Jo {s — v) 


l-u 


dv 


= D^-^Zg. 


Therefore, 


Consequently, 


Yt = 


(t) + ^ D'^nt - s)D^-’^Zgds 


= Yf 


Yt = ^ r{t - u)y^gdBl. 


A Technical appendix 

In this section, we gather some useful results from [38] and recall a theorem on the regularity 
of the convolution product. We denote by F[^ the set of functions on [0,1] with Holder 
regularity A. 
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A.l Fractional integrals and derivatives 

Lemma 13.1 in [38] relates the Holder exponent of a function and the Holder exponent of its 
fractional derivatives. 

Proposition A.l. If f £ and /(O) = 0, then for any a < X, f admits a fractional 
derivative of order a and D^f £ 

Equation 2.20 in [38] is a fractional integration by parts formula which can be written as 
follows. 


Proposition A.2. If cf £ and with 1/p 1/q < 1 + a, then (f and ip have an 

integral of order a and 


[ (j){t — s)I^'ip{s)ds = [ I°‘(l){t 
Jo Jo 


s)^p{s)ds. 


In this work, we mainly use the two following corollaries of Proposition A.2. 


Corollary A.l. Let (p £ , with r > 1 and ip £ . Then, for any a < /3, D°‘ip exists, 

belongs to and 


f (p{t — s)ip{s)ds = f I'^(p{t — s)D°‘ip{s)ds. 
Jo Jo 


Corollary A.2. Let (p be continuous and ip such that x^ip{x) £ for some p > 0. Then, 
for any a < min(l — p,X), D°‘ip exists, belongs to L^ for some r > 1 and 


[ (p{t — s)ip{s)ds = f I°‘cp{t — s)D°‘ip(s)ds. 

Jo Jo 


A. 2 Convolution 

Finally, the next result is about the smoothness of the convolution of a power type function 
with a continuous function. 


Proposition A.3. Let f be a differentiable function on (0,1] such that for some K > Q, 
0 < /3 < 1 and any x in (0,1], 

and g a continuous function on [0,1] . Then the convolution 

f *g{t) = [ f{t- s)g{s)ds 
Jo 


has Holder regularity (1 — /3). 


18 



Proof. We write G for the supremum of \g\ and we split / * g{t + h) — f * g{t) into the three 
following terms: 


f*git + h)-f* g{t) = 


ct-\-h 


f{t + h- s)g{s)ds 


+ [ + s) - f{t-s))g{s)ds 

J t—h 

+ [ {f{t + h-s)-f{t-s))g{s)ds. 

Jo 


The first term is bounded by the second by KG{1 + j^)h^ ^ and the third by 


pt—h pt-\-h—s p 

/ / \f{u)\duds <GK 

Jo Jt—s Jo 


G 




^ ds < ^GKh^-^. 


□ 
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